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The short- and long-time dynamics of model systems undergoing a glass transition with apparent 
inversion of Kauzmann and dynamical arrest glass transition lines is investigated. These models 
belong to the class of the spherical mean-field approximation of a spin-1 model with p-body quenched 
disordered interaction, with p > 2, termed spherical Blume-Emery-Griffiths models. Depending on 
temperature and chemical potential the system is found in a paramagnetic or in a glassy phase 
and the transition between these phases can be of a different nature. In specific regions of the 
phase diagram coexistence of low density and high density paramagnets can occur, as well as the 
coexistence of spin-glass and paramagnetic phases. The exact static solution for the glassy phase is 
known to be obtained by the one-step replica symmetry breaking ansatz. Different scenarios arise 
for both the dynamic and the thermodynamic transitions. These include: (i) the usual random 
first- order transition (Kauzmann-like) for mean-held glasses preceded by a dynamic transition, (ii) 
a thermodynamic first-order transition with phase coexistence and latent heat and (iii) a regime of 
apparent inversion of static transition line and dynamic transition lines, the latter defined as a non¬ 
zero complexity line. The latter inversion, though, turns out to be preceded by a novel dynamical 
arrest line at higher temperature. Crossover between different regimes is analyzed by solving mode 
coupling theory equations throughout the space of external thermodynamic parameters and the 
relationship with the underlying statics is discussed. 


I. INTRODUCTION 

In the present work we investigate the dynamic prop¬ 
erties of a glassy system in which, under certain exter¬ 
nal conditions, both glass and fluid can coexists, yielding 
different scenarios for dynamical arrest and for the fluid- 
glass transition. These properties can be studied in sta¬ 
tistical mechanical models with bosonic spin-1 variables, 
where the holes s = 0 play the role of inactive states, that 
is, the so-called Blume-Capel [1, 2] or Blume-Emery Grif¬ 
fiths (BEG) [3] models. In these models the fluid phase 
corresponds to a paramagnet and the solid phase is ei¬ 
ther a ferromagnet (no or weak disorder) [1-5] or a spin 
glass (strong disorder) [6-12], In the present work we 
consider an extension to p-spin interacting systems with 
spin-1, to p > 2 and continuous (spherical) variables [13] 
to better represent continuous density fluctuations, alike 
to liquid-like compounds. 

In the presence of quenched disorder the random BEG 
model with pairwise (p = 2), as well as its spheri¬ 
cal counterpart, is known to display both a continuous 
paramagnet/spin-glass phase transition and a first-order 
one (first order in the thermodynamic sense, i.e. with 
latent heat and a region of phase coexistence). Further¬ 
more, melting upon cooling [14-20] can occur, with a spin 
glass at high T and a paramagnet at low T. These prop¬ 
erties have been observed in the mean-field approxima¬ 
tion, where the self-consistent solution for the spin-glass 
phase is computed in the full replica symmetry breaking 


(RSB) Parisi ansatz [21] and on the cubic 3D lattice with 
nearest-neighbor couplings [22, 23]. The frustrated BEG 
model has been studied by means of numerical renormal¬ 
ization group techniques, as well, with results depending 
on the underlying lattice and the renormalization tech¬ 
nique adopted [24-26]. 

Mean-Held spin-glass models with Ising [27], soft [28, 
29] or spherical [30, 31] spins with more than two-spin 
interactions, called p-spin models, are known to yield 
the so-called random first-order transition, i.e., a phase 
transition across which no internal energy discontinuity 
occurs but the order parameter (the Edwards-Anderson 
overlap Qea) jumps from zero to a finite value. Their 
glassy phase is described by an ansatz with one RSB 
[32]. In a cooling procedure, the thermodynamic tran¬ 
sition is preceded by a dynamic transition due to the 
onset of a very large number of metastable states sepa¬ 
rated by high barriers [33]. “Very large” means that the 
number of states A f grows exponentially with the size N 
of the system: AT ~ exp(ElV) where the coefficient E is 
the configurational entropy, also called complexity in the 
framework of spin-glass systems (see, e.g, Refs. [34, 35] 
and references therein). “High barriers” means that the 
free energy difference between a local minimum in the 
free energy functional of the configurational space (also 
called free energy landscape) and a nearby maximum (or 
saddle) grows with N. The phenomenology of the p-spin 
spin-glass systems is, in many respects, very similar to 
the one of structural glasses. These models are, there¬ 
fore, sometimes called mean-field glasses. The occurrence 
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of non-zero £ is a fundamental property both in mean- 
field systems [28, 29, 36] and outside the range of valid¬ 
ity of mean-field theory, e.g. in computer glass models 
[37, 38], or, indirectly, by measuring the excess entropy 
of glasses in experiments, see, e.g., Ref. [33] and refer¬ 
ences therein. The barriers’ height turns out to diverge 
in the thermodynamic limit in the mean-field approxi¬ 
mation, this being an artifact of mean-field glasses. The 
thus induced dynamic transition corresponds to the tran¬ 
sition predicted by another mean-field theory for the dy¬ 
namics of supercooled liquids: the mode coupling theory 
[39]. The thermodynamic transition occurring at a lower 
temperature is, instead, the mean-field equivalent of the 
so-called Kauzmann transition in glasses, also known as 
the ideal glass transition [33]. This was initially predicted 
by Gibbs and Di Marzio [40] and its occurrence in real 
strutural glasses is still object of an ongoing debate [41 
44], 

We are going to investigate the complex dynamic prop¬ 
erties consequent to the combination of a kind of interac¬ 
tion inducing structural glass behavior and the presence 
of hole states (aka, spin state s = 0) inducing phase co¬ 
existence. The latter element is, possibly, responsible for 
melting upon cooling [21, 22]. The first of such mod¬ 
els was brought about by Sellitto in the pairwise ran¬ 
dom orthogonal model with spin-1 variables [45]. In the 
present dynamic work we rather consider the multi-body 
interaction model of Ref. [13], where both high temper¬ 
ature coexistence of high- and low-density paramagnetic 
phases, and low temperature coexistence of (low-density) 
paramagnetic and spin-glass phases are displayed. 


II. MODEL 

The model we consider is a spherical Blume-Capel [1, 2] 
model with p-body disordered interactions. Our starting 
point is the model Hamiltonian 

H = — ^ ■ ■ ■ s ip T D 'y ] Si , (1) 

ii<i2<---<ip i 


lattice-gas 

ii<i2<---<ip 

+ (D — Tlog 2) ^ iii, 

i 

where the —Tlog2 term is necessary to keep the ratio of 
filled-in to empty sites identical to the one of the original 
Hamiltonian, see [47, 48] for details. We then introduce 
the variable [49] 

Tj = o-j(2n, - 1) = ±1. (3) 

This way, the model Hamiltonian assumes the form 

H ~ — 2P 'y 1 Ji 1 i2-..ip{ cr ii + Til) . . ■ ( a ip + T i P ) 

ii<i2<---<i P 

+ (T-Tlog2)^^^, 

i =1 

( 4 ) 

where all the degrees of freedom, cr’s and r’s are now Ising 
spins. A continuous spin model can then be constructed 
by imposing two independent spherical constraints 

N N 

= J2 t ' = n - ( 5 ) 

i-1 i=1 

Programma/Pspinbosonic/Test/ The thermodynamic 
properties of the model we just defined are thoroughly 
studied, applying 1RSB theory, in reference [13]. For the 
sake of brevity, henceforth we shall refer to Ref. [13] as 
the static study (that is, replica theory-based), in con¬ 
trast with the dynamic study that constitutes the prin¬ 
cipal aim and subject of this paper. Before deriving and 
solving the dynamical equations of the model, it is suit¬ 
able to briefly summarize the static results of [13], with 
particular emphasis on the aspects that will be most rele¬ 
vant for the dynamical study that we are going to report. 

A. The static phase diagram 


where the variables s, are bosonic spins (i.e., Sj = 
1,0,—1), and the couplings are independent quenched 
random variables distributed with a Gaussian probabil¬ 
ity density, 


P(Ji 1 i 2 ...i p ) — 


i Np- 


nj 2 p\ 


exp 


I 2 

J lll2- 


NP~ 


p\ 


The 1RSB free energy for the model is 


B 2 

f3F =(/3D — log2)d —— (d p + ( m — l)9i) — log 2 

1 f Tfl — 1 1 

- O 1 °g( 1 ~ d ) + —-logr/o 4-!°g »7i 

2 V m m 


( 6 ) 


the external parameter D is called “crystal field” in lit¬ 
erature, and it essentially plays the role of a chemical 
potential. Because of the bosonic spins, we cannot de¬ 
fine the continous spin approximation of the model with 
spherical constraint in the usual way [46]. We must first 
rewrite the Hamiltonian as an Ising-spin problem on a 


with the definitions 

Vo = d — q\, (7a) 

Vi = d+ (m - l)qi - mq 0 , (7b) 

where the parameter d is the ratio of filled-in to empty 
sites (that is, the density of the system), and the q 0 and 
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gi are respectively the mutual overlap and self-overlap, 
as usual in a 1RSB Ansatz. The extremization of the (6) 
whith respect to gi and d yields the saddle-point equa¬ 
tions 


P/3 2 

2 


dP~ l 


P/3 2 p-i 

‘r q i 

pP 2 p -i 

-2-rf 


gl 

pom ’ 


(8a) 


Vo — @ 
Oi 7o 


+ 2(/3£>-log2). (8b) 


The static lines of the phase diagram, in absence of 
an external magnetic field, can be determined by setting 
go = 0 and studying the system given by the equations 
(8a), (8b), and the saddle-point condition for the 1RSB 
parameter m 


z(y) = 


( 9 ) 


where z(y) is the Crisanti-Sommers function[30] 


z(y) = -2 y 


1 - y + log y 

(1 - y) 2 ’ 


_ Vo 

y = — 

Vi 


= A (10) 


The dynamical arrest lines can be, as well, identified us¬ 
ing replica theory, by solving the saddle-point equations 
for m = 1, and by imposing the following marginality 
condition for the solution in the d , gi parameter space: 


d 2 F 

a 2 f 

dq{ 

dqidd 

d 2 F 

a 2 F 

dddqi 

dd 2 


This corresponds to looking for a spinodal point in the 
free energy as a function of d and gi. 

The resulting phase diagram for the model is reported 
in figure 1. 

From it, it can be seen that the system exhibits a 
rich phenomenology, with both random first-order transi¬ 
tions (RFOT) and thermodynamic first-order phase tran¬ 
sitions (TFOPT). Here we will just comment on them 
briefly to adequately introduce our dynamical study, 
along with statics novelties with respect to previous anal¬ 
ysis. The interested reader can find all the details in [13]. 

a. Random first-order transition For low enough D, 
the system exhibits the random first-order phenomenol¬ 
ogy typical of the p-spin model. Along the dynamic tran¬ 
sition line Td(D), the system undergoes a dynamical ar¬ 
rest , meaning that the relaxation towards the paramag¬ 
netic, stable state is blocked by the presence of an ex¬ 
ponentially large number of metastable SG states which 
trap the dynamics. Being the SG states metastable, this 
transition is not captured by the static saddle-point equa¬ 
tions, and has to be studied by solving the dynamics of 
the model, or by using the marginality condition for the 
dynamics [30, 47, 50, 51]. 

At a temperature T S {D) lower than Td(D), a static 
transition takes place, whereupon the number of SG 
states becomes subexponential (their complexity, also 



FIG. 1. The complete phase diagram for the model (4) for 
p = 3. In the bottom inset the detail contained in the box of 
the main panel is displayed, together with he values of Da, 
D t h e -Dinv (see text) identified by vertical dotted lines. 


called configurational entropy, vanishes) and they be¬ 
come stable with respect to the paramagnet, yielding 
the equilibrium SG solution [30, 52]. This feature of the 
model, occurring for D < Db = 0.77, is equivalent to the 
p-spin model phenomenology. 

b. Thermodynamic first-order phase transitions In 
the region of the phase diagram between the spin¬ 
odal lines, but above the RFOT line, two paramagnetic 
phases, termed PM - and PM + , coexist, both with gi = 0 
but with two different density values d: d~ and d + > d~. 
The two paramagnets are labelled + and — according to 
their density value d being, respectively, large and small. 
These values can be determined by solving the saddle- 
point Eq. (8b) in the gi = 0 limit, yielding the expres¬ 
sion 

^d p (l — d) = T 2 (2d— l) + 2(TD — T 2 log2)(l —d)d. (12) 

For p = 3, this is a polynomial equation with three 
solutions for d £ [0 : 1]. The solution with the inter¬ 
mediate value of d turns out to be always unstable (see 
Ref. [13]), leaving only a high-density and a low-density 
solutions. Since the density d has a continuous behav¬ 
ior along the PM—/SG thermodynamic transition, Eq. 
(12) can be used, as well, to determine the value of the 
density for the SG phase at the transition point. This 
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means that the spinodal lines T(D) can be determined of the Hamiltonian (4) from the deterministic one 
in a parametric form in d as 


D sp (d) = 


d S 


£-1 


(13) 


S pv~v Y p - 1 2(d - 1)VM 2 - Ad + 2 

|d 3 (p - 1)2 log(2) + d 2 [p(2 - 41og(2)) + 41og(2)] 
+d [p(log(4) - 3) + 1 - 2 log(2)] + p| 


H = H 0 + Hj, 

N , i 

H 0 = (D-Tlo e 2)<T i ™^, 

Hj ~~~ 2 p . E Jiit2 ••• ip(Oil + Til) ... ( ^~ip)' 

ii<i2< mmm <ip 


Tsp(d) 


{d- 1 )\/p(p- 1 )d p/2 
V4d 2 - 4d + 2 


(14) 


The relaxation dynamics is, then, governed by the 2N 
Langevin equations 


By studying those expressions, it can be readily checked 
that for p > 3, at high D the spinodal curve is an asymp¬ 
tote of the D axis; this means that, however large D , the 
system can always present an high density phase, if the 
temperature is low enough. This fact will be of capital 
importance in the following of this paper. 

Besides the RFOT, thus, the system also exhibits both 
a PM—/PM+ and, furthermore, a PM—/SG thermody¬ 
namic first order transition, that means standard first 
order transitions with phase coexistence and latent heat. 
Both transitions take place along the TFOPT line in fig¬ 
ure 1 and coexisting phases exist between the spinodal 
lines. As the temperature is decreased to cross the RFOT 
dynamic line, we observe that only the high-density para- 
magnet undergoes dynamical arrest, while the low den¬ 
sity PM~ phase is unperturbed. Thus, at T lower than 
the crossing point of TFOPT and RFOT lines, the tran¬ 
sition occurs between a low density paramagnet PM—, 
with d = <i_ and qi = 0, and an lrigh-density SG with 
d = d + and q± ^ 0. The intersection takes place for 
(D b ,T b ) = (0.77,0.40). 

c. High density dynamical transition. At 
(Z} t h, T t h) ~ (0.876,0.379), the dynamical RFOT 

line and the spinodal SG line intersect. It can then 
be seen that for D > D t h, the dynamical RFOT line 
coincides with the spinodal line of the TFOPT, which 
means that the dynamical arrest in the PM + phase will 
take place as soon as phase separation occurs. From 
the thermodynamic point of view, we have coexistence 
between two paramagnets, as before. However, if we 
perform a quenching dynamics from the high density 
PM phase, a dynamical arrest into a metastable SG 
phase will take place. 


III. THE DYNAMICS 


We are now ready to derive the dynamical equations 
for the model. Let us first separate the disordered part 


&i = -n{t)<Ji{t ) 
Ti = -v(t)Ti(t) 


dH 0 

d<Ji 

dH 0 

dn 


dHj 

d<7i 


+ Vi(t), 


dHj 

dn 


+ 9i(t), 


(15) 


where we assume the noise fields rji(t) and 9i(t) to be 
delta-correlated: 


= 6ijD 0 (t - t'), 
{Vi(t)0j(t')) = 0, 


with (taking the Boltzmann constant ks = 1) 


Do(t — t') =2 T6(t-t'). 


Following [53], we have inserted the Lagrange multipliers 
and i/(t)Ti(t ) in order to enforce the spherical 
constraint. 

The quantities we are interested in are the correlation 
functions and the response functions of the system. In 
this case, differently from the standard spherical p-spin 
model, we have two different types of degrees of freedom, 
and thus two possible external perturbing fields, one for 
each of them. As a result of this, we have four different 
correlation functions 


Caa(t,t’) = -^^(cr i (t)cr J (t , )) ) 

i=1 

1 N _ 

Cari^f) = — ^ (o-i(t)Ti(t')), 

V i =1 

i N _ 

Cra(t,t') = 

i =1 

1 N _ 

Crr{t,t') = — ^ (n(t)Ti(t 0 ), 
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and four response functions 


here the site-independent dynamical equations 


7? (t t' s 1 1 V 

’ ’ N Z—i ShAt') ’ 

a = —n<j(t) 


+K p (p 

R (t 1 V' 

i—1 v 7 

+£(*)> 

Rr,(t,t)- , 

x=l y ' 

Ti = —vr(t) 

R it 1 V 

Rrr(t,t) - , 

+K p (p 

i—1 

+C {t)i 


{D-T log 2) 


it) 


(16a) 


(D-T log 2) 


it) 


(16b) 


where hi(t) and li(t) are time-dependent perturbing 
helds conjugated with the a, and r, degrees of freedom, 
respectively. 

Our aim is to use equations (15) to obtain self- 
consistent dynamical equations for the functions above. 
We employ the generating functional method devised in 
[54] by Martin, Siggia and Rose, and already used for the 
p-spin model in [55] by Kirkpatrick and Thirumalai. The 
MSR approach consist essentially in defining a generating 
functional Z[h, l, h, l, J^] [56] for the 2N 1-dimensional 
random fields tJj(t) and iy(t)[57]. The correlation and re¬ 
sponse functions can, then, be obtained by taking func¬ 
tional derivatives of Z with respect to the external fields 
h(t), l(t ), h(t) and l(t), as in an usual field theory. 

We have emphasized the fact that the generating func¬ 
tional still depends on the quenched random couplings 
Jij , and so does every quantity generated by it; so, in 
principle, we would have to average them over the disor¬ 
der in order to obtain the correlation and response func¬ 
tions we want. However, as remarked by De Dominicis 
in [58], since the generating functional in absence of ex¬ 
ternal currents is by definition normalized to one 


where the correlation matrix for the noise terms £(f) and 
Q(t) has been renormalized in the following way 

(Zi(t)Sjtf)) =6 ij K p C(t 1 ty- 1 + 6ij2TS(t - t')), 

= «i(t)<jtf)) (17a) 

UtMt')) =S ij K p C(t,t l r~ 1 . 

We have also defined the constant 
K =^ 

1X P ~ 2 2p+l 

and the two functions 

C = C a(7 + C aT + CV CT + C' TT , (18a) 

7Z = Raa + Rot + Rto T Rtt- (18b) 

The equations for the correlation and response function 
can then be easily obtained, as reported in appendix A. 

A. Symmetries, equilibrium and ergodicity 


Z[0,0,0,0,Jy] = 1, 

it is independent from the variables of the system, and so 
it can be averaged over the disorder directly. This is in 
contrast with the static partition function for a system 
with quenched disorder, which is not self-averaging. As 
a result of this, the use of replica theory is not needed in 
the dynamical framework; this fact constitutes the main 
advantage of the dynamical approach over the static one. 

Performing the average over the disorder leads to a de¬ 
coupling of the lattice sites and a coupling of the configu¬ 
rations of the system at different times, as it happens for 
the p-spin model in [55]. This is to be conceptually com¬ 
pared to the decoupling of the sites in the static approach 
with replicas, yielding a coupling between different repli¬ 
cas [59]. It is then possible, using saddle point methods, 
to write an effective generating functional which yields 
two single site dynamic equations valid for every degree 
of freedom in the lattice. The details of the derivation of 
the dynamics can be found in the appendix. We report 


In appendix A, we derive eight coupled differential 
equations for eight different unknown functions. We now 
specify them to the particular problem we want to study, 
i.e. identifying dynamical arrest. In order to to this, 
we can restrict ourselves to an equilibrium (i.e. starting 
from an equilibrium initial condition) and ergoclic dy¬ 
namics. This implies time-translational invariance (TTI) 
of the correlators 

C(t, t!) = C(t — t') (19) 

and that the fluctuation-dissipation theorem (FDT) 
holds 

R(t) = -9(t)^^-, ( 20 ) 

where R and C denote any correlation and response func¬ 
tion couple in the system, respectively, and 9(t) is the 
Heaviside step function. These assumptions are valid in 
the high temperature PM phase, where ergodicity is not 
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broken, but they are generally false when the system is 
cooled below the dynamical transition temperature Td, 
where a transition to a SG phase with broken ergodic- 
ity takes place. Second, we notice that both the model 
Hamiltonian (4) and the effective generating functional 
(A9) are symmetric with respect to a c f> t switch 

H(cr —► r, r —» <x) = H (<x, r). 

This means that the er and t evolve in the same statisti¬ 
cal ensemble, which, in turn, implies that the correlation 
functions obey the relations 

c aa (t,t') = c TT (t,t'), c aT (t,t') = c Ta {ut'), (21) 

that can then be extended to the response functions by 
exploiting the FDT 

R„„ = = Rrr, ( 22 ) 

R °* = = Rt °- (23) 


which is a closed equation, whose solution can be found 
once the values of the parameters D , T and d are known. 
It can be seen as a Mode Coupling schematic equation, 
and it is worth noticing that in the limit of the d param¬ 
eter (the fraction of filled-in sites) going to 1 one recovers 
standard MC equations. The density d is assumed, at all 
times, to be equal to its equilibrium value, given by Eq. 
(12). As a result of this, in the region between spinodal 
lines, where phase coexistence occurs, the dynamics will 
have to be studied separately for each one of the two co¬ 
existing phases. In the high density limit case d —> 1, 
thermodynamically occurring for low D , the proper limit 
of Eq. (27) is considered in appendix A3. The study of 
Eq. (27) will be one of the main focuses of this work. 

An equation can be derived also for the difference AC 
between C aa and C aT , obtaining 


d(AC(t)) 

dt 


(D-T\og2)d-^dP-T\ AC(t), (28) 


which has the trivial solution 


Once that these relations are established, we can see 
that only the two correlation functions C aa and C aT are 
needed to completely describe the dynamics of the sys¬ 
tem. Thus, we can define the “total” correlation function 
of the system C(t — t'), by normalizing the C to one 


A C{t) = 2(1 - d)e At (29) 

where the A constant is defined as 

A = (D — T log 2)d — -j^dP — T. 


^ C(i-0 2 ( Caa (t-t') + C ar (t~t')) 

[ 1 C(0) 4d(0) 

where we have used the spherical constraint 


N 

E (Ci(t)o-i(f)} = IV => C a<T (t,t) = 1 

i=1 


and the relation (3) between the err product and the site 
occupation number 


Vi{t)Ti(t) + 1 
2 


rii(t) 


C aT (t,t) 


C TtT (t,t) = 2 d(t) - 1. 


(24) 

(25) 


Since we assume the dynamics to be at equilibrium, the 
density d(t), being an one-time quantity, is a constant of 
motion, equal to its equilibrium value 


d(0) = d(t) = d( oo) Vi. (26) 


Using the equations for the C aa and C aT , we obtain 
the following integro-differential equation for the total 
correlation function C(t) 


C(t) = A(D,T,d)C(t) 

—B(T, d) f du C p ~ 1 {t - 
Jo 

{ 1 2 HP~ 1 T) 

A(D,T,d)= \ D-Tlog2-—^- 


(27) 

u)C(u) 

) (d-l)-T 


B(T,d) 


J 2 pdP~ l 

4T 


If now we use Eq. (12) to eliminate D, we get 


2(1 - d) 

which is always negative for any value of d and T. This 
means that the difference between the two correlators 
tends to zero for long enough times, and thus the dy¬ 
namics of the system can be always solved using the C(t ) 
function only. 


B. Solving the dynamic equations 

In this section we report the results obtained by nu¬ 
merically solving equation (27) in various representative 
points of the phase diagram. Eq. (27) is an integro- 
differential mode-coupling like equation [39], that can be 
solved using the standard algorithm intruced by Fuchs et 
al. [60] and extended in different ways, cf. e. g., Refs. 
[61-63]. We consider three meaningful cases to illustrate 
the varoius occuring regimes. 


1. D = 0: dynamic transition 

At D = 0 and high T the system yields a single param¬ 
agnetic phase (PM), which undergoes a dynamical tran¬ 
sition with ergodicity breaking at Td{D = 0) = 0.4892 
[13]. In figure 2 we report the total correlation function 
of the system as the transition is approached from above. 
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FIG. 2. The total correlation function C{t) evaluated for var¬ 
ious temperatures in the vicinity of the dynamical transition 
at Td(D = 0) = 0.4892; for each temperature, the value of the 
equilibrium density d of the system is reported. As the sys¬ 
tem is cooled, the function develops a plateau whose length 
increases rapidly with decreasing T, to become infinite at the 
dynamical transition line. 


FIG. 3. The total correlation function C(t) evaluated for var¬ 
ious temperatures as the system enters the phase coexistence 
zone and the high density phase undergoes dynamical arrest. 
Here Td(D = 0.77) ~ 0.40586. In the phase coexistence zone, 
the PM+ correlator is plotted with a thin line, while the PM— 
with a thick one. As we can see, the PM+ makes a transition 
to an arrested SG phase while the PM— remains paramag¬ 
netic. 


As we can see, our dynamic equation (27) yields the 
dynamic transition predicted in [13], and the correlator 
shows the typical phenomenology of a mode-coupling like 
dynamic arrest. The transition temperature Td corre¬ 
sponds up to a O(10 -4 ) error with the predicted value 
0.4892 and the height of the plateau q,i (also called the 
non-ergodicity parameter) is (p — 2)/(p — 1) = 0.5, as 
expected. 


2. D = Db: phase coexistence 

For D = Db = 0.77, the situation is richer and more 
interesting; for this value of the crystal field, at low tem¬ 
perature the system undergoes phase coexistence, yield¬ 
ing two separate paramagnetic phases with high (PM + ) 
and low (PM - ) density. As we anticipated, the dynamics 
of the system has to be solved separately for each one of 
these two phases. Their behavior turns out, actually, to 
be quite different, as only the high density PM+ phase 
undergoes ergodicity breaking as the dynamical line is 
crossed. In figure 3 we plot the resulting correlators as 
the system enters the phase coexistence zone, and the 
high density phase undergoes the dynamic transition. 

We see that, again, the expected phenomenology is 
reproduced by our equation. At T = 0.42586 only a sin¬ 
gle paramagnetic phase is present, but for T = 0.41586, 
two different paramagnetic phases separate; at T = 
0.40586 ~ Td(D = 0.77), the high density phase correla¬ 
tor shows the typical plateau as the dynamic transition is 
approached, and as the cooling continues, the high den¬ 
sity PM + phase eventually undergoes a dynamical arrest, 


while the low density paramagnet PM - stays unchanged 
until the first order phase transition takes place to the 
thermodynamic IRSB-stable glass phase. 


3. D > Dth-' anomalous dynamical arrest 

Up to the prevuoius case, our dynamic analysis con¬ 
firms the results of the static one [13]. However, as we 
anticipated in section I, for D > D t h , our dynamical 
equation yield dynamic arrest at finite temperature when 
starting from the PM + phase, an effect that was not iden¬ 
tified in the static analysis, where the thermodynamically 
dominant phase is low density PM - . We shall now report 
the results obtained by solving equation (27), leaving the 
study and discussion of the new arrested phase for the 
next section. 

In order to show the onset of this new transition, we 
choose D = 1 > D th , T = 0.35 and we cool the system 
starting from a high density initial condition crossing the 
spinodal line during the procedure. The PM + spinodal 
line, for these chemical potential values, turns to be a 
dynamical arrest transition line, as shown in figure 4. We 
can observe that, as soon as the two PM phases separate, 
the high density phase is already an arrested SG phase, 
with a nonzero overlap, while the low density phase shows 
no sign of dynamical arrest. This is at difference with 
respect to the static results of Ref. [13], where the high 
density phase is supposed to be still paramagnetic for 
T = 0.34. 

What is puzzling about this result is the fact that the 













high density phase is already deep into the SG when the 
separation occurs: for T = 0.34 we have qd ~ 0.82, which 
is already much higher than the (p — 2)/{p— 1) = 0.5 that 
we would expect for a system which approaches from 
above an usual dynamical arrest. This means that if 
the high density phase existed even above the spinodal 
line, its dynamic transition temperature would be actu¬ 
ally much higher than T = 0.34; however, this effect is 
not visible since only the low density PM~ phase exists 
in that region. If we follow the same cooling procedure 
for D = 2 (i.e., we cross the spinodal line at D = 2 dur¬ 
ing the cooling), the results are not very different, and 
the high density arrested phase is still present. In this 
case, the value of the overlap at the separation line is 
even higher, with q,i ps 0.97 for T = 0.2. 

Solving equation (27) for higher D does not change 
the general situation, so we will not report any results 
for higher values. The point is that, since the spinodal 
line is an asymptote of the D axis (as we have men¬ 
tioned in section II), then for arbitrarily large D , a high 
density phase exists at T > 0. According to both the 
marginality condition for the statics, and the dynamic 
results presented so far, this phase presents a dynamical 
arrest into a SG phase with nonzero overlap, that can be 
realized by selecting atypically dense initial conditions at 
those temperature and chemical potential. 


IV. COMPLEXITY AND FREE ENERGY 


Since the system undergoes a dynamical arrest, we 
would expect the RFOT phenomenology which holds in 
the other regions of the phase diagram to be present in 
this case as well, inside the high d minimum which corre¬ 
sponds to the PM+ phase (the low density paramagnet 
is completely orthogonal to our discussion). In summary, 
we expect the metastable states which trap the dynamics 
(and maximize the complexity) to have a higher in-state 
free-energy than the one of the paramagnetic, ergodic 
state. We might also expect a complexity £(/) to be 
strictly positive for every / up to a static temperature 
T a {D) where states with null complexity are born and a 
static transition takes place. We computed both quan¬ 
tities using replica theory as in Ref. [13] and report the 
corresponding curves in figures 5 and 6. We see that 
the usual RFOT picture does not hold for all values of 
D along the new dynamical arrest line. Indeed, it holds 
only for D < Di nv where Di nv ~ 0.9062, which in fact 
corresponds to the point where the static line touches 
the new dynamical arrest line. In this region, the com¬ 
plexity of the states is positive for every / (as one can 
see in the upper panel of figure 6) and the free energy of 
the m = 1 states which trap the dynamics is higher than 
the paramagnetic one. On lowering T at constant D, the 
complexity touches zero (corresponding to the condition 
Eq. (9)) and a static transition to a SG phase takes place 
as usual. 

This picture does not hold anymore for D > Di nv — 
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FIG. 4. The total correlation function C(t) evaluated for var¬ 
ious temperatures as the system enters the phase coexistence 
zone and the anomalous dynamic transition takes place. The 
presence of this dynamic transition is a novelty with respect 
to the results of [13]. 


0.9062. For D = Di nv , the maximum complexity con¬ 
dition (dynamical arrest line) and the zero complexity 
condition (static line) coincide, which means that the 
complexity is zero (i.e., the number of all states be¬ 
comes subexponential) and the free energy of the trap¬ 
ping states becomes equal to the paramagnetic (PM+), 
cf. figure 5. This is alike to the occurrence of the static 
Kauzmann transition but in this case both SG and PM + 
phases are metastable and dynamically occurring only 
for initial conditions with a density atypically high for 
these chemical potential D values. When D is increased 
along the high density arrest line, the free energy of the 
trapping states becomes lower than the paramagnetic one 
and continues to decrease indefinitely, while the complex¬ 
ity becomes more and more negative as reported in fig¬ 
ure 6. If one looks at the phase diagram 1, it is possible 
to see that there is a range of D such that, on cool¬ 
ing, the static line is met before the dynamic one [13]. 
This would mean that along this path in the parame- 
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FIG. 5. The state free energy of the low density paramagnet 
(PM-), the high density paramagnet (PM+), and of the spin 
glass states which trap the dynamics (SG), evaluated along 
the high density dynamical arrest line for D > D t .h = 0.876. 
We also report the densities of the two paramagnets. We see 
that the free energy / of the trapping states becomes equal 
to the one of the paramagnet for D = Di nv ~ 0.9062 


ter space the states with low free energy and complexity 
arise before the ones with maximum complexity and free 
energy, only eventually appearing on the maximum com¬ 
plexity line (that, for D < D t h, is the dynamic line): the 
complexity of the glassy metastable states becomes again 
non-zero crossing the static line on cooling and reaches 
its maximum at the line formerly denoted as dynamic 
line. However, this is true for D > D t h — 0.876 up until 
D ~ 1.6. From that point on, the static line touches the 
T = 0 axis and the complexity stays negative for every 
value of D,T. In summary, using both the dynamical 
equations and the replica approach, we find an anoma¬ 
lous, complexity free, dynamical transition, that occurs 
if the system is initially prepared at values of the den¬ 
sity corresponding to a paramagnetic metastable state 
at high density. This state always coexists with a more 
probable and thermodynamically dominant low density 
paramagnetic state. It is the latter which, at lower T 
and/or lower D , undergoes a first order phase transition 
to a spin-glass phase. 



f 



FIG. 6. The complexity of the states which trap the dynamics, 
below the high density dynamical arrest line. Upper Panel: 
complexity for two points with D < Di nv ~ 0.9062, with 
T S (D) < T < T d (D) = ... (full line) and T = T S {D) = ... 
(dotted line), respectively. Lower panel: complexity for two 
representative points with D > Di nv . 


V. CONCLUSIONS 

In this work, we have studied, both dynamically and 
statically, a disordered model which shows both an 
RFOT-like phenomenology (dynamical arrest, complex¬ 
ity, Kauzmann transition, etc.), and an “ordinary” first 
order phase transition with latent heat and phase coex¬ 
istence. We have derived and solved the equations for 
the equilibrium dynamics of the model and completed 
the static replica-based study of [13] with novel results. 
In doing so, we have noticed the presence, in a certain 
region of the phase diagram, of a RFOT-like dynamical 
arrest line, which however shows a non-positive complex¬ 
ity and does not work as a precursor for a Kauzmann-like 
static transition, in contrast with the usual phenomenol¬ 
ogy expected in RFOT models. 

The picture we propose to explain this is the following. 
Our model has two order parameters, the density d and 
the 1RSB self-overlap q± (termed q from now on). One 
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can imagine to construct a potential function of d and q , 
by plotting the paramagnetic free energy as a function of 
d , and then performing a Franz-Parisi like [50] construc¬ 
tion along the q axis, for every d. For D < Db = 0.77 
the transition scenario is qualitatively alike to the ran¬ 
dom first order transition one. 

For chemical potential values in the interval Db < 
D < D t h , instead, the scenario changes. For high T only 
one minimum is present but, on cooling, the spinodal 
line of a second paramagnetic phase, at higher density, 
is crossed: a second, metastable minimum with q = 0 
and density d = d + is formed, termed PM + . Further 
lowering the temperature the dynamical transition line 
is crossed. There, a metastable phase with g / 0 and 
d = d + arises from the PM + solution, with a higher free 
energy. This phase consists of many equivalent states 
and corresponds to the arrested glassy phase. The differ¬ 
ence in free energy between the SG metastable minima 
and the free energy of the PM + minimum is equal to 
the complexity counting the log of the number of SG 
metastable states. If T is lowered further, we cross the 
static transition line, where the q ^ 0 minimum has the 
same height as the paramagnetic one (null complexity) 
and thus it becomes stable: we have a static transition 
in a spin-glass phase. 

To summarize, we have a three-step process on cooling 
for D e [D b , Ah]: 

1. A secondary PM + minimum with d = d+ and q = 0 
is formed (the stable phase is the low density PM - ). 

2. A SG minimum with d = d+, q ^ 0 is formed, 
arising from PM + . 

3. The SG becomes stable with respect to the low den¬ 
sity paramagnet PM - and a static transition takes 
place. 

This scenario almost corresponds the usual RFOT phe¬ 
nomenology. The only difference is that there are 
two paramagnetic phases, one stable (PM - ) and one 
metastable at higher density, PM + , and that the glassy 
metastable states at the threshold free energy arise inside 
the PM+. 

When D is raised beyond Ah) the order of these steps 
above is changed. Two apart scenarios appear. For 
Ah < D < Anv) steps 1 and 2 exchange. Fixing T 
slightly above T sp (D), d = d + and looking at the poten¬ 
tial only in the direction of q , the g / 0 minimum will 
have already formed. However, this has no effect on the 
thermodynamics since d = d + does not yet correspond 
to a minimum on the d axis. For that to happen, T must 
be lowered to cross the PM+ spinodal. At lower T the 
process goes on as previously. 

For D > Am/, step 1 becomes the last to happen: 
again, the system dynamics is arrested as soon as the 
PM + spinodal is reached, but the SG minimum along q 
not only is formed before the PM + has a chance to ap¬ 
pear, but it has even become stable with respect to the 
PM + . We stress the fact that this unusual behavior is 


possible only because the system has two order parame¬ 
ters, differently from usual RFOT models whose behavior 
is governed only by q. The fact that the minimum in q is 
already formed for D > _D t h is, indeed, confirmed by the 
fact that the naive marginality condition d 2 F/dq 2 = 0, 
used in Ref. [13], completely misses this new arrest line: 
in the direction of q the minimum is already formed and 
the curvature is positive. The whole Hessian matrix of 
the replicated free energy with respect to q and d (eval¬ 
uated at m = 1) must be used to detect the new line, as 
reported in Sec. II A. This picture is corroborated by the 
fact that this phenomenology is found in the PM + phase, 
whose density goes up and approaches 1 as D increases, 
as reported in the upper panel of figure 5: this means 
that the PM + phase of our model becomes more and 
more similar to the usual p-spin spherical model (PSM) 
[30], and the dynamical equations behave accordingly as 
explained in appendix A3. However, we also stress that 
such density values are thermodynamically extremely un¬ 
likely to occur for these large values of the chemical po¬ 
tential D. 
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Appendix A: Derivation of the dynamical equations 
for the p-spin Blume-Capel spherical model 


1. Computation of the effective generating 
functional 


In this appendix we provide the step-by-step derivation 
of the dynamical equations studied in this paper. 


We write down the full expression for the MSR gener¬ 
ating functional 


Z[h, l, h , f, Jij] = / Va'Da'Dr'Dr x exp 


— - / dtdt!cr(t) ■ <r(t)D 0 (t — t') + f(t) ■ T(t')D 0 (t — t') 


x exp 


x exp 


x exp 


J dt i ( ff(f) • (<x (f) + p(t)cr(t) + d a H 0 (t )) + f • (f(f) + v{t)r{t) + d T H 0 (t ))) 

J dt i (<r(f) • dcrHj(t)) + f (t) • d T Hj(t)) 

— f dt i(h(t ) • d-(t) + • cr(t) + Z(t) • r(f) + l(t) ■ r(t)) 


(Al) 


where V denotes a functional integration measure over disorder is contained only in the Hj part of the liamilto- 
all lattice sites; for example nian, we have to perform the average only on the third 

line of the (Al). Thus we have to compute the integral 

N 

Va = 6cr(t). 

i—l 

The first step is to perform the average over the disorder 
of expression (Al) for the generating functional; since the 


exp 


J dt i (&(t) • dcrHjit)) + f(t) • d T Hj(t)) 

f dJi 1 ir 


exp • 


J z . . N p ~ l t r 

"l! l 2 ...lp- L ' °l\l2--^p / 7 1 rv- - - \/ , \ / | \ 

/ dt [lycr^ + + r i2 )... (a ip + n ) 


J 2 p \ 2 2 Pp\ 

\ 

+*(cr n + r il )( l j i2 + f i2 ).. .(<j ip + Tj p ) H-h {a tl + T h ){a i2 + n 2 )... i(d ip + T ip )\ 


(A2) 


where we have symmetrized the (<r + f)^ (<r + r)i 2 ... (a + r)i couplings; performing the gaussian integration yields 


exp 


J dt i (<x(f) • d a Hj(t)) + f (f) • d T Hj(t)) 

J dtdt'^p i(& + t )fi(cr + f) t > (cr + r) t • (tr + ■ 


22p+2jYp-l 


p-l 

t' 


r i r i r i p_ ^i 

p{p - 1) i(& + r)t ■ (a-+ r) t ' (<t + r) t -i{cr + r) t > (tr + r) t ■ (cr + r) t > >■ 

(A3) 


As we anticipated in section III, the average over the 
disorder has decoupled the lattice sites, at the price of 
generating a coupling between configurations of the sys¬ 
tem at different times. 


This dynamical coupling is conceptually similar to the 
coupling between replicas that occurs in the static treat¬ 
ment of the p-spin model [30, 52]; so, following [55], we 
define the dynamical overlaps between auxiliary MSR 
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fields, 


N 


1 

i=l 

Qg(t,t') 

1 N 

= j^y 

l— 1 

Qio(t,t') 

1 w 

Qz{t,t') = —yiTi(t)iai(f), 

i=l 

Qn(t, t') 

1 w 

QaM = — yiTi(t)iTi(t'), 

Ql2{t, t') 


and between auxiliary and dynamical fields 
1 N 

Qd(t, t') = — y t) 

V i= 1 
1 N 

Qlo(M') = Jy ^2iVi(t)Ti(t') =Qii(t',t) 

i=1 
1 N 

Qn(t,t') = jyy^ifj(f)(7i(f') =Qi 5 (t',t) 

i =1 

1 * 

Ql2(t,t') = —^iTi{t)Ti{t') =Ql 6 (t',t) 


i=l 


between dynamical fields, 


AT 




These functions can be incorporated in the (A3) by in¬ 
serting 1-factors in the form of functional Dirac deltas 
expressed in exponential form, for example 


*=l 
N 


Qe{t,t') = — <Ji(t)Ti(t'), 


1 = / SQiSXi exp 


i =1 
N 


QrW) = 


dtdt'iXi(t, t') 
x | NQi(t,t') - y i<7k(t)iak{t') 


N 


i=l 

N 


(A4) 


<3s(i, t ) — ^ ' T i{t)Ti{t ), 


the (A3) then assumes the form 


J S 16 QS m Xexpl^J dtdt'N y iX/jQ/, 
x exp <| 22 p +2 J dtdt' p(Qi + Qi + Qz + + Qq + Q 7 + 

p(p — 1) ^9 + QlO + Qll + Ql2^ ^Ql3 + QlA + Ql5 + Ql6^ + ^6 + Ql + 

x exp i J dtdt' [X\ia ■ ia + Xgicr ■ if H-+ X\qT • if] | . 


p-2-i 


This expression contains an exponent proportional to the 
system size N, which in the thermodynamic limit allows 
us to perforin the functional integrals over the Qs and As 
using the saddle-point method. The saddle-point equa¬ 


tions read 


iXi... 2 A 4 — — 


tA 5 ...iAs = — 


J 2 p 

22p+2 

J 2 p{p - 1) 


(Qs + Qe + ^ + W 1 , (A5) 

(Qi + Q'2 + Qz + Qa) (A6) 


2 2 p+ 2 

x (Q 5 + Qe + Q 7 + Qs) p 

J 2 p{p- l)(p- 2) 


2 2 P+2 

x (Ql3 + Ql4 + Q 15 + Qw) 
(Q.5 + Qe + Qr + Qs) p 3 ) 


(Q 9 + Qio + Q 11 + Q 12 ) 
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J^viv — I") 

*Ag . . . iX\2 = -22p+2-(Q 13 < 5 14 + Q 15 + Qw) 

x(Q 5 + Q 6 + Q7 + Qs) p 2 1 (A 7 ) 

7)(t) — 

*Ai2 ■ ■ • *Ai6 =- 2 2 P+ 2 -^ 10 

x (Q5 + Q 6 + Q7 + < 58 ) p 2 - (A 8 ) 

The Qi -4 can be self-consistently set to zero, as they 
essentially are the correlation functions for the auxiliary 
MSR fields; besides this, it is immediate to see that the 
second term in the (A6) contains the sum of all the re¬ 
sponse functions for the system Qg_i 2 , multiplied for a 
sum of the same functions with inverted times Q 12 - 16 ; 
thus, since every response function R(t,t') is zero for 
t < t' because of causality, the resulting product van¬ 
ishes and the Ag_g are all equal to zero. 

Using the saddle-point equations, we can eliminate the 
As and replace the Qs in the saddle point with the corre¬ 
lation functions, 

= c aa (t,t'), 

Qi p (t,f) = c aT (t,t'), 

Q^ p (t,t') = C Ta (t,t'), 

Ql P (t,t') = C TT (t,t'), 


and the response functions of the system 

Qf p (t,t') = -R aa (t,t , ) = Ql P (t , ,t) 

Ql4 P (M') = Raripi t') = Qfo (t',t), 

Qf[ (tj) = -Rra{t,if) = Qf P (t',t) 
Qf[(t,t') = -RrrM = Q S 12 P (t',t). 


If we now define 


K = _— 

p ~ 2 2 P + 1 ’ 

and 


C — C aa + C aT + C Tr7 + C TT , 

/v = R a u + R GT ~\~ R TG H - RtT-) 

we can rewrite the averaged generating functional in the 
following way 


Z[ 0] = J VgVcfVtVt exp 
x exp J dt id 
x exp [ dt if 


rr + per 


J dtdt'(dD 0 d + fD 0 f) exp — ^ J dtdt'K p (d + f)C(t, t') p 1 (d + f 

(D-T\og2) T _ _ 1} f dt , t , )c ^ t y-2 {a{tl) + T(t / )} 


(A9) 


{D-T log 2) 


■ — K p (p— 1) J dt' TZ(t,t')C(t,t') p 2 (cr(t') + r(t')) 


as in [55], this functional is now local in space (i.e. refers 
to single sites only) but nonlocal in time; we can also see 
that this expression is quite similar to the one we had in 
equation (Al): we have the quadratic terms in d and r, 
containing the noise correlators, and the linear ones con¬ 
taining the equations of motion (16a), (16b) themselves. 


2. Equations for the correlation and the response 


Once the effective generating functional has been com¬ 
puted, the self-consistency equations for the correlation 
and response functions can be readily derived. For the 
Cuer we just use the basic definition 


dCuu(ti,t 2 ) 

dti 


(d(fi)o-(f 2 )), 


where (•) denotes the average over the gaussian thermal 
noises £(£) and £(f)[64]; using the (16a), we get 

dc **(h,h) = _^ tl)C(j(j{ t ut2) _ OP-rlog2 ) c T u(t 1 ,t 2 ) 


dti 


+K p (p- 1 ) J dt'n(t,t')C(t,t') p - 2 

x {C a u(t ', t 2 ) + C T u(t', t 2 )) + (£(ti)cr(t 2 )). 


The last term has to be computed using directly the 
general definition of an MSR generating functional (see 
[52] for details). First we write 

Z = z N 


where the z is the reduced generating functional for a 
single couple of dynamic variables (a , r) only. With this 
definition we can write (£(ti)cr(t 2 )) as (to lighten the no- 
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tation, we omit the time contractions). 


for the C a<7 


fafafafa ) = H)tt 


5£ 5( 5a 5cr 5 t 5t 


r 

x cr fa) exp 
x exp [i(— £<t + £7 - < 0)1 exp [£] 


3=0. 

-\(ZG €€ Z + 2ZGkC + CG k Q 


(A10) 


where the probability distribution of the noise has been 
made explicit by defining the inverse of the correlation 
matrix 


G = D~\ 


D K {tfa) = mcfa)), 

D iC (t,t') = m«t')). 


dC cra fa 1 fa 

dt\ 


llfa)C acr fa,fa 
_ (D-Tlog2 )_ CrgihM) 

fa 

+ K p (p -1) dt' lZ(t, t')C(t, t’f - 2 

J — 00 

x (Caafa 1 ^ 2 ) + C Tr7 fa, fa) 

+ 2TR aa fa, fa + K p r dt'Cfat'y- 1 

J —OO 

x (R (7(7 t') + -R £7T (^2 5 0); 

(A13) 


the term proportional to the is zero since we as¬ 
sume t\ > t 2 - The same method can be applied to the 
remaining correlation functions. 

For what concerns response functions, it can be easily 
proven that[52, 65] 


5afa)\ _ 5 (a fa)) 
5^ fa) / 5hfa) 


(A14) 


and all the terms that don’t contain the noise fields (the 
time-derivatives of the dynamic variables, the Lagrange 
multipliers, the convolutions with the response and cor¬ 
relation functions) have been cropped in C. 

The functional integral over the noise is a standard Gaus¬ 
sian integral with a linear term, so it can be readily per¬ 
formed, yielding 


J 5a5d5r5T a fa) exp j - ^ [(-<r + j)Dfa-a + j) 

+ 2(-<t +j)Dtf(-T) + tD cc t] | exp [C] ; 

(All) 


now, by taking the functional derivative with respect to 
jfa) and multiplying by —i we get 


{fafaafa)) = J 5a5d5r5T J dt! [D^(ti, t')(—i)afa)afa) 
+ Dfa tl ,t')(-i)afa)ffa)]e s ^^ 

= J dt'[D^fa,t')R a<7 fa, t') 

+ D^fa,t')R (7 T fa,t')\, 

(A12) 


where e' s '( cr - T ’ ,T - T ) is the probability distribution for the dy¬ 
namic and auxiliary variables induced by the probability 
distribution for the noise [54]. 

By using this result, and recalling the definition (17a) for 
the correlation matrix, we can finally write the equation 


which in turn implies 


dR aa fa,fa _ Jfa I Safa) \ _ / Safa) \ 

dti dt 1 \5^fa)/ \5€(t 2 )/’ 

and so 

dR aa fa,fa / 5a{t i)\ 

(D-T log 2) / 5rfa) \ 

2 \5^fa)/ 

fa 

+ K p (p -1) / dt 1 nfa t')c{t, t'y~ 2 

J — 00 


Safa) 

Mfa) 


Srfa) \ 

Kfa)/_ 


+ 


m 1 ) 
sfafa 

(A15) 


The other equations can be obtained by trivially gener¬ 
alizing the (A14). 

Let us make an important remark. The last term in the 
(A15) (given by the functional derivative of the noise with 
respect to itself) is a Dirac delta, which assures that 


lim R a <jfa,fa = I- 

1 1 — 2 " 


The same applies to R TT . On the other hand, the equa¬ 
tion for the R aT is the following 


dRgrfafa) 

dti 


— nfa)R aT fa,fa — - - - - ^-^-Rcrcrfafa) 

fa 

+ K p (p -1) dt' Ufa t')C(t, t'y “ 2 

J — 00 


(R 

<JT fa , fa + R (7(7 fa, fa), 


(A16) 
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which indeed lacks the Dirac delta, implying 


lim R aT (t 1 ,t 2 ) = 0, 

t x —^2" 


which applies to the R Ta too. This relation is funda¬ 
mental to recover the static equation for the density (see 
appendix A 4). 

The last thing left is computing the expression for the 
Lagrange multipliers p(t ) and z/(t); since we introduced 


those terms to enforce the spherical constraints, their ex¬ 
pression must be determined self-consistently from them, 


C acr {t , t) = N 
C TT (t,t) = N => 


dC aa (t,t) 

dt 

dC TT (t,t) 

dt 


= 0 , 
= 0. 


By expressing these total derivatives in terms of the 
dynamical equations for the correlation functions, we get, 
for the p(t ) 


and for the v(t), 


V(t) 


+ K P (p- 1) f dt'n{t,t')C(t,t') p - 2 {C^(t',t) + C T(7 (t',t)) 

J — OO 

+ TR tt (£, t) + Kp f dt'C(t,t') P 1 {Raa (£> t') H - R(TT (t-) t')) — 0, 
J — OO 


( t) = _ ( D-rio S 2) (2 „ w _ 1) 

+ K p (p - 1) [ dt'K(t, t')C{t , t'r 2 (C TT (t', t) + C aT {t ', <)) 

+ TR Ga (t, t) + Ap f dt'C(t,t') p 1 (Rtt{ t, t') + R Tr7 {t, t 1 )) = 0. 

J—00 


(A17) 


(A18) 


These are the expressions for /u(f) and that are to be 
plugged in the equations for the correlation and response 
functions. 

From this point, the derivation proceeds by imposing 
the FDT to couple the response and correlation func¬ 
tions, and then using the TTI to manipulate the resulting 
expressions. Those steps are only a matter of standard 
algebra, so we shall not report them. 


3. The d —> 1 limit 

The Blume Capel p-spin spherical model is supposed 
to be indistinguishable for a plain PSM in the d —> 1 
limit, since in this limit the sites are all filled-in. So, our 
equations are supposed to give back the classic dynamical 
equation for the PSM when d —> 1. 

From equation (28), we can see that C<uff{t) is auto¬ 
matically zero for all times when d=l. So we only need 
to take care of equation (27), but the problem here is 
more tricky. In fact, it is immediate to verify that by 
just sending d to one in the (27), we get a wrong result. 
The problem is that the value of d depends strongly on 
the point (D,T) of the phase plane by means of equation 
(12); indeed, if we invert the picture and we treat d as a 
free parameter, we see that when d —> 1, then the abso¬ 


lute value of the crystal held D tends to infinity. This is 
perfectly reasonable; when D becomes lower and lower, 
the configurations with a good number of empty sites are 
energetically disadvantaged, and the partition function is 
dominated mainly by configurations with many hlled-in 
sites, which become the only ones accessible to the sys¬ 
tem when D — 00 . Since \D\ —> 00 at the same time 
as d approaches one, the latter limit is nontrivial and has 
to be taken carefully. 

Therefore, in order to estimate the order of magnitude of 
D as d tends to one, we solve the (12) for D , obtaining 

3 d p+1 - 3 d p + d 2 T 2 log 16 - 4T 2 d(log 2 - 1) - 2 T 2 
~ 4(d- l)Td ‘ 


By plugging this into equation (27), we have 


C(t) = 


T 

T - T 

2d 


J 2 p(d)P~ 1 r* 

4 T In 


C(t) 

du C ,p_1 (t — u)C(u); 


by imposing d = 1, this becomes 

C(t) = ~C(t) - ^ f du C p ~\t - u)C(u), 
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which is the dynamical equation for the PSM [52, 53, 55], 
apart from a factor 1/2 in front of the r.h.s.. Its presence 
is due to the fact that in the limit d —> 1, a = r and 
so our model contains 2 N a spins, while the variance of 
the J couplings has the usual scaling for a system of N 
spins. By taking the long time limit, we get 


g = p P 2 a r>~ 1 

1 - q 2 q 


(A19) 


which yields the well-established values for the dynamic 
transition temperature T d and the dynamical overlap q d 


Qd 


P~ 2 
P~ 1 


T d 


p(p — 2 )p~ 2 

2 (p- 1)p- 1- 


(A20) 


This shows that, since in the long time limit C(t) is 
always zero, the 1/2 factor is completely irrelevant for 
what concerns the properties of the system at the tran¬ 
sition point. 


4. The equation for the density 

In this work, we have stressed multiple times the fact 
that the static equation (12) is needed to solve the dy¬ 
namics of the system, since the equilibrium value of d 
is not known a priori. This is fine for all practical pur¬ 
poses, but clearly unsatisfying from the theoretical point 
of view, as the (12) has been derived using replica theory 
in a static scenario. It is then necessary to show that 
equation (12) can be derived also in the MSR formalism 
in order to get a completely auto-consistent theoretical 
picture. 

Our starting point is a trivial rephrasing of relation (3) 

n, Vi(t)Ti(t) + 1 

=-2-' 

We sum this relation over lattice sites, divide by the sys¬ 
tem size N and then take the thermal average on both 
sides. We get the following relation between the density 
and the err correlation function 

«*» = (A21) 

By taking the total derivative with respect to time we 
get 


equations for the C aT 


d(d(t)) 

dt 


= -M(*)( 2 (<*(*)> - 1) - 


D-T log 2 


J' 2 p{p - 1) 


dt 1 1z(t, t'y 


'\p- 2 


2 2 p+l 

x (C<7T (t', t) + C TT {t !, t )) + 2 TR ra (t, t) 


(A23) 


Now, given that lim ti _^+ Ra T {ti,t 2 ) = 0, and using the 
TTI and FDT, we get 


d(d(t)} 

dt 


p(t)(2d(t)-l)^(^ D ^° g2 ) 

T 2 n 2 r°° 

2^r J Q c ^y~ lc ^ du ' 


(A24) 


where we have switched to the notation (d{t)) = d(t). 
The expression (A17) for the Lagrange multiplier can be 
manipulated in the same way, obtaining 


d(t) = - D T 2 l ° g2 (2d(t)-l) 
k r°° 

" 2 # J du C(u)C(u) p ~ 1 +T; 


(A25) 


note that, since we are looking for a dynamical equation, 
we are now working with a time-dependent d. 

We can now see that the integrals in both those expres¬ 
sions can be easily computed using the substitution for¬ 
mula. By doing so and than plugging the expression for 
p in the equation we get 


d(t) 


^W)f(2rfW-i) 

+ ^ f mt)r-T(2d(t)-i), 

(A26) 


which, after some trivial algebra, leads to 

d = 2 d(d - 1)(D -T log 2) - JJpd p (d - 1) - T{2d - 1). 

(A27) 

If we now assume the system to be at equilibrium at all 
times 


d{d(t)) 1 dC(t,t) 
dt 2 dt 

= 1 f dC ( rr(t 1 ,t 2 ) dC <7T (t 1 ,t 2 )\ 

2 ^ dh dt 2 

ti=t 2 =t 

(A22) 

which can be then rewritten in terms of the dynamical 


d = 0 d(t) = d eq 

the (A27) can be rewritten as 

|d p ( 1 - d) = T 2 (2d - 1) + 2 (TD - T 2 log 2)d(l - d), 

which indeed is the static equation (12) for the density. 
This is the result we wanted. 
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At this point, the reader has probably noticed the fact 
that, in this derivation, we have used the TTI and FDT, 
which are suited to an equilibrium situation, to derive a 
dynamical equation for the thermal average of a one-time 
observable, which is supposed to vary only in a nonequi¬ 
librium scenario. 

This indeed is not a problem in itself. We could have em¬ 
braced the equilibrium picture from the very beginning 
(i.e., setting dC ^' t ' > = 0 in the first place), and we would 
have anyway recovered the static equation for the den¬ 
sity, which was the only result we needed for our dynam¬ 
ical study to be self-contained. However, the dynamic 
equation (A27) can be interesting to study in itself; by 
completing it with an initial value for the density d(0), 
we obtain a Cauchy problem for which an unique solu¬ 
tion is guaranteed to exist. We plot four such solutions 
in figures 7 and 8 for (D,T) =(2,0.2). 


(D,T) =(2,0.2) 

d(t) 



(D,T)=(2,0.2) 

d(t) 



FIG. 7. Solutions of the (A27) for ( D,T ) = (0.2,2), with 
d(0) = 0.1 (up), e d(0) = 0.8 (down); both solutions converge 
rapidly to their respective equilibrium values, d” = 0.051100 
and d+ = 0.934117. 


The plots show that, in the phase coexistence region, 
both equilibrium solutions correspond to an attractor for 
the dynamics, which is not surprising; where the solution 


actually ends up depends on the initial value d(0) that is 
chosen. 

Besides this, it is interesting to notice that the unstable 
root dmt of equation (12) (which, as we said in section 
II, is always to be discarded on grounds of replica sta¬ 
bility arguments) is not an attractor for the dynamics, 
and actually is the critical value of d(0) that marks the 
separation point between the two attractors. Thus, the 
selection rule of section II for the value of the density is 
recovered from a dynamical point of view. 


(D,T)=(2,0.2) 

d(t) 



(D,T)=(2,0.2) 

d(t) 



FIG. 8. Solutions of the (A27) for (D,T) = (0.2,2), with 
d(0) = 0.754865. slightly below dint (up), and d(0) = 
0.754867 slightly above di„t (down), dint = 0.754866. The so¬ 
lution ends up in a different attractor depending on the choice 
of the initial value d(0). 


We conclude this section with a final remark. It may 
seem contradictory that the density d(t) evolves in time, 
since it is a one-time quantity which is assumed to be 
constant in an equilibrium dynamics scenario. However, 
this is true only if the initial value d(0) (which is com¬ 
pletely arbitrary) is chosen different from both of the two 
equilibrium solutions of the static equation, which is ex¬ 
actly the contrary of what condition (26) requires. 

Thus, the (26) can be seen both as a reasonable physical 
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assumption, and as a self-consistency condition for the equilibrium dynamics. 
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